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ABSTRACT 

We present the results of microlens ray-tracing simulations showing the effect of ab- 
sorbing material between a source quasar and a lensing galaxy in a gravitational lens 
system. We find that, in addition to brightness fluctuations due to microlensing, the 
strength of the absorption line relative to the continuum varies with time, with the 
properties of the variations depending on the structure of the absorbing material. We 
conclude that such variations will be measurable via UV spectroscopy of image A of 
the gravitationally lensed quasar Q2237+0305 if the Lyman Alpha clouds between the 
quasar and the lensing galaxy possess structure on scales smaller than ~ 0.1 pc. The 
time scale for the variations is on the order of order years to decades, although very 
short term variability can occur. While the Lyman alpha lines may not be accessible 
at all wavelengths, this approach is applicable to any absorption system, including 
metal lines. 

Key words: gravitational lensing - quasars: individual: Q2237+0305 - quasars: ab- 
sorption lines - intergalactic medium 



1 INTRODUCTION 

The Lyman Alpha (Ly Q ) forest is an absorption phe- 
nomenon which is observed in the spectra of many high 
redshift quasars. It arises due to absorption by galactic and 
intergalac tic neutral hy drogen between a quasar and an ob- 
server [see lRauchl il99 8'l for a recent review] . As the universe 
expands, photons from the quasar are redshifted, and there- 
fore absorption by clouds at different cosmological distances 
cause absorption at different wavelengths in the spectrum. 
The resulting spectrum contains many absorption lines blue- 
ward of the Ly a emission line of the quasar, and the struc- 
ture of the absorbing materia l along the line of sight is im- 
printed on the spectrum fe.g. lLu et al.lll99r?) . 

Numerical simulations of cosmological structure for- 
mation predict the large scale structure of the gas which 
can then be compar ed with observations of quasar spec- 
tra jDave et al.lll999h . The conclusions of these studies are 
that the Ly a "clouds" are mostly intergalactic structures, 
with the denser regions corresponding to the locations where 
galaxies form. The absorbing gas along the line of sight to 
a single quasar can be probed by comparing the relative 
positions and strengths of the absorp tion lines in the Ly a 
forest in the spectrum of the quasar ijSargent et alJ ll980L 
The structure across the line of sight is harder to measure, 
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although if two quasars that are nearby on the sky are ob- 
served, the fraction of Ly a forest lines that are common to 
the two quasars can pr ovide a measure of the extent of the 
clouds across the s ky jSargent. Young, fc Schneider! 11983 : 
iRollinde et afll2003lb 

Traditionally, the large scale structure of the Ly a For- 
est has been the main focu s of research, because of its 
cosm ological significance (e.g. lMcDonald fc Miralda-Escudl 
1999). Recently, the small scale structure (on kilopar- 



sec scales) has also been investi; 


;ated by several methods 


( Rauch. Sargent, fc Barlow! [1999 , 


200l| iRauch et all 


2001, 


2002); potentially revealing structure in and about the in- 



terstellar medium in young galaxies. 

Gravitational lensing has previously been used to probe 
the size of absorbing clouds [note that if the mass asso- 
ciated with a Ly a cloud is sufficiently l arge, it too ma y 
result in gravitational lensing effects, see iMenardl <2004h ]. 
When a background quasar (regarded as a point source) 
is lensed by a foreground galaxy, multiple images of the 
quasar may be obtained. The light rays responsible for the 
two different images have traversed slightly different regions 
of space in their journey to the observer. Hence, a measure 
of the size of the Ly a absorbing clouds (down to kiloparsec 
scales) can be made by cross-correlating the the absorption 
lines in two images of the same quasar fellison et alJl2004l : 
iBecker. Sargent, fc Rauchll2003ri . This method is analogous 
to using a true pair of quasars, however the lensing supplies 
a "pair" of quasars that are close together on the sky (sepa- 
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Figure 1. As the relative positions of the observer, lens and 
quasar change (eg due to rotation and linear motion of the 
galaxy), the ray passes through a different region within the ab- 
sorption clouds; causing the absorption line strength to change 
with time. This is s imilar to the changing positions of the micro- 
imag es with time iLewis fc Ibatal Il998t iTrever fc WambsgansJ 

rated b y ~ arc seconds), allowing smaller scales to be inves- 
tigated. |Rau£li J _SM^ejit 1 _&_£arlow ( 1999J) used this method 
in the multiply imaged quasar Q1422+231 to trace Mgll 
and CIV absorption clouds between the lens and the quasar 
down to scales of ~ 26 pc, which is the smallest scale so far 
considered. 

The very small scale (< 1 parsec) structure of the Ly a 
clouds is completely unknown (due to limits in observational 
techniques, and the fact that the resolution of Ly a forest 
formation simulations is too coarse), and would be useful 
for deciding what the local equivalents of the Ly a clouds 
are. Some (but not all) low redshift L y a clouds have been 
found to be associated with galaxies (Salp eter fc Hoffman! 
Il995l) . suggesting that the structure of the Ly a ckmds might 
mimic that of the interstellar medium. This is known to have 
structure on many scales, exhib iting fractal structure due 
to turbulence felmegreenl I199 1 ??) . In this paper we present 
simulations demonstrating how the very small, sub parsec 
scales of Ly Q clouds can be measured using a quasar that is 
subject to gravitational microlensing by a foreground galaxy. 



2 BACKGROUND AND APPROACH 

In gravitational microlensing, the lensing effects of in- 
divid ual star s in th e lensing galaxy are considered [see 
IWambseanssI il998l) for a review of gravitational lens- 
ing, including microlensing of quasars]. The cluster of 
' micro-images' that result are too small to be resolved 
jLewis fc Ibatall998tlTrever fc Wambsganssl2004l) . with the 
quasar appearing point-like, however brightness fluctuations 
can result as the relative positions of the source, observer 
and lensing stars change; this is the observational conse- 
quence of microlensing. Since the rays joining the observer 
to the quasar traverse different regions of space at different 
times, we anticipate that small scale structure in the absorp- 
tion clouds will lead to time variations in the strength of the 
absorption line relative to the continuum (see Figure 0. 

Throughout this paper, we assume a cosmology with 
Q m = 0.3, Ov = 0.7 and H = 72 km s _1 Mpc -1 . 

2.1 Method 

We m ade use of the microlensing code of IWambseanssI 
(1990), modified to output the positions of the rays in the 
lens plane and the corresponding position in the source 
plane. The code uses the backward ray tracing method, 



where a rectangular grid of rays is fired into the lens plane, 
and the corresponding position in the source plane is cal- 
culated by summing the microlensing effects of many stars 
in the lensing galaxy. Using these positions, magnification 
maps (plots of the magnification of a point source as a func- 
tion of position in the source plane) can be produced by 
binning the positions of the rays in the source plane (ie 
areas which collect a lot of rays have high magnification). 
For the simulations presented in this paper, typically > 10 9 
rays were traced between the observer and the source, en- 
suring that the resulting magnification maps, consisting of 
1024x1024 binned pixels, did not suffer significant noise from 
Poisson sampling. A third plane, the "cloud plane" was in- 
troduced, positioned between the lens plane and the source 
plane. 

For any ray starting at a position (a;, y) in the lens plane 
and landing at a position (x a ,y s ) in the source plane, the 
position (x p , y p ) at which that ray intersects the cloud plane 
is simply found by linear interpolation: 

x p = x + X(x s - x) . . 

Vv = v + Hys -y) 

where A is a parameter which determines the location 
of the cloud plane. If A = 0, the cloud plane coincides with 
the lens plane, and can be use d to model absorption clouds 
within the lensing galaxy (e.g. iLewis fc Ibatall2003^l . Alter- 
natively, A = 1 corresponds to the case where the absorption 
clouds are located at the source. This situation has been used 
to investigate the effects of microlens ing on broad absorp- 
tion line quasars iLewis fc Belldll99a) . Note that 
are scaled distances, but can be interpreted as the angular 
position of a point measured in units of the angular Einstein 
radius 

_ / 4GMq D la 

e °-y^^D^D7 s (2) 

Hence, if a point in the cloud plane has scaled coordinates 
{x p ,y p ), the physical coordinates (eg. measured in parsecs) 
are (D op 9ox p , D pQoy p ), where D op is the angular diameter 
distance from the observer to the cloud plane. 

The clouds were described by an absorption function 
A(x p ,y p ) defined over the cloud plane, which is the frac- 
tion of photons that are absorbed as a function of where 
the ray intersects the cloud plane. If a ray passes through 
a region where A = 0, it is unaffected, but if a ray passes 
through a region where A = 0.5 say, then the ray is atten- 
uated by 50%. Gaussian clouds of the form A = \eT T ^ b 
were used, distributed randomly over the cloud plane. To 
investigate the effect of the cloud size b, we ran simulations 
with three different cloud sizes for each value of A. The cloud 
sizes used were b — 0.005 pc, b — 0.02 pc and b = 0.1 pc 
(for typical cloud configurations, see Figure|2j. The number 
of clouds in the cloud plane was chosen so that the aver- 
age absorption per unit (physical) area of the cloud plane 
was approximately the same for all of the simulations (this 
was not exact due to a 'clipping' effect, since the absorption 
can never be greater than 1 where several clouds coincide 
and their absorption values add) . The mean value of the ab- 
sorption was ~ 0.5. To achieve computational efficiency, the 
absorption as a function of cloud plane position was stored 
as a 1024x1024 pixel image, and A(x p , y p ) was evaluated for 
each ray by looking up the appropriate pixel of the image. 
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Figure 2. The three different absorption cloud sizes which were 
investigated. The cloud sizes are 0.005 pc, 0.02 pc and 0.1 pc 
respectively. The plots show the absorption as a function of posi- 
tion in the cloud plane. Black corresponds to complete absorption 
(A = 1), white corresponds to no absorption (A = 0). Note that 
due to the external shear, the rays actually pass through a rect- 
angular area in these images (sec Figure These images of the 
clouds were for A = 0.2, the physical size of the plane was different 
in the other two cases to ensure that all the rays passed through 
the cloud plane. 



The physical size of the cloud plane was chosen to be large 
enough to ensure that all of the rays intersect it and can be 
affected by absorption. 

For the purposes of this study, the microlens- 
ing parameters adopted were those of image A of 
Q2237+0305, with a dimensionless surface mass density 
(all in stars of m ass 1 Mr) of a = 0.36 and exter- 
nal shear 7 = 0.41 jWa mbsgans s. Paczvnski. fc KatzllT990i 



ISchmidt. Webster, fc LewislQ998i) . The region of the source 
plane that was considered was a square with side length 20 
Einstein radii, ie a physical length 209oD os , however this 
ray bundle is only square at the source plane. It is ap- 
proximately rectangular for the rest of the journey, with 
side lengths as a function of redshift shown in Figure 
which was calculated assuming that the lensing galaxy was 
composed of smoothly distributed matter (see section l2~2t . 
The redshift of the lens is z — 0.0394 and the redshift 
of the source is z = 1.695, which are the redsh i fts fo r 
Q2237+0305 and its lensing galaxy dHuchra et allll98Eft . 



For these redshifts, the physical size of an Einstein radius 
in the source plane is 5.85xl0 -2 pc and the speed of the 
quasar across the source plane was assumed to be 600 

km s -1 , or approximately one Einstein radius per d ecade 
jKavser. Refsdal. fc Stabe3ll986t iLewis fc Ibatalli99sh . We 
ignored the individual motions of stars in the lensing galaxy 
which, for a transverse velocity of 600 km s _1 would cause a 
10 per cent increase in the frequency of high magnification 
events llKundic fc WambsganssUl993T) . Ignoring this effect 
saves considerable computational effort. 

The magnification maps provide a measure of the inten- 
sity of the source as a function of its position in the source 
plane. However, for a non-point source, the total observed 
intensity is the source profile multiplied by the magnification 
map, integrated over the source plane. Therefore the inten- 
sity as a function of source plane position for an extended 
source is given by the convolution of the magnification map 
with the source profile. The physical size of the continuum 
sou rce in Q2237+0305 has been measured to be < 7.8xl0 -4 
pc iShalvapin et al.l 120021) . which is ~ 0.25 pixels wide in 
the magnification maps. Therefore, we ignored finite source 
effects, taking a single pixel to represent a point-like source. 

For each simulation, a pair of magnification maps was 
produced. One was the usual magnification map calculated 
using the final locations of the rays, neglecting absorption, 
representing the magnification of the unabsorbed (contin- 
uum) wavelengths of the quasar spectrum. An absorbed 
magnification map was also produced, using the attenuated 
values of the rays. This magnification map applies to the 
absorbed part of the spectrum, so measures the flux at the 
bottom of the absorption line as a function of source po- 
sition. A third map was produced, a "line strength map", 
which measures the flux at the bottom of the continuum 
relative to the flux of the continuum as a function of source 
position. It is obtained by dividing the absorbed magnifica- 
tion map by the unabsorbed map (see Figure 0. 

A sample light curve showing the flux in the continuum 
and the absorption line, as well as the strength of the absorp- 
tion line as functions of time is shown in Figure|l| This light 
curve was taken from the same magnification maps shown in 
Figure |3 by taking a one dimensional horizontal cut across 
the middle of the magnification map. Note that some (but 
not all) of the sudden changes in continuum brightness (due 
to the source crossing a caustic in the source plane) are ac- 
companied by sudden changes in the line strength. Also, for 
times when the continuum magnification is approximately 
constant (eg from times ~ 0.6 — 1.4 in Figure 0J, the line 
strength time series may still vary, reflecting the structure 
in the absorption clouds. 

2.2 Distances and Cosmology 

The scaled lens mapping equations (which relate a position 
in the lens plane, or the image, to a position in the source 
plane) are 

x a = x — atx(x,y) 
y s = y~a y (x,y) 



(3) 



where a x and a y are the scaled deflection angles. Substi- 
tuting these equations into Equation Q converting back to 
unsealed coordinate s and using the unsealed lens equation 
JWambsgansslll99d) . it can be shown that A is given by 
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Figure 3. Examples of magnification maps showing the magnifi- 
cation as a function of position of the source (the cloud plane po- 
sition was A = 0.5 and the cloud size was 0.02 pc for these plots). 
The first is the unabsorbed magnification map for the continuum, 
the second is the absorbed magnification map for the flux in the 
hy a absorption line. The third is the ratio of the two, and mea- 
sures the strength of the absorption line as a function of quasar 
position (light areas denote a weak absorption line, dark areas 
indicate a deep absorption line). 
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Figure 4. A typical light curve showing how the flux in the ab- 
sorption line (dashed line) and the continuum (solid line) vary 
with time (produced by taking a horizontal track across the mag- 
nification maps). The structure of the absorption clouds produces 
variations in the line flux relative to the continuum (dotted line). 
The time units depend on the velocity of the quasar, for a typical 
velocity of a few hundred km s~ x this plot spans a time scale of 
several decades for Q2237+0305. The y-axis was normalised to a 
mean of 1, corresponding to the theoretical mean magnification, 

Vth = ( 1 _ (T )2_ 7 2 • 
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Figure 5. The interpolation parameter A as a function of the 
redshift of the cloud plane for Q2237+0305 (* ; = 0.0394, z s = 
1.695). 



where D oa , Di v , D op and Di s are the angular diameter dis- 
tances from the observer to source plane, lens plane to cloud 
plane, observer to cloud plane and lens plane to source plane 
respectively. The redshifts were converted to angular diam- 
eter distances using the assumed cosmology. The resulting 
plot of A vs z is shown in figure |5] Note that the flatness of 
this relationship over a wide redshift range implies that the 
results for A = 0.95 are approximately applicable for clouds 
over a very broad redshift range. Hence, the very small scales 
of the Ly a forest could be probed over a wide redshift range 
(z ~ 0.5 to z ~ 1.6) by comparing observations with the 
high A simulations. Also, the physical size X p of any feature 
in the cloud plane with scaled size x p is simply 

X p = D op 9 x p (5) 

In this paper, we investigate the cases A = 0.2, A = 0.5 and 
A = 0.95. For these values of A, the corresponding redshifts 
for the clouds are z = 0.049, 0.077 and 0.53 (although the 
latter result would also work well for clouds at higher red- 
shifts). Absorption by neutral hydrogen at these redshifts 
would produce an absorption line in the Q2237+0305 spec- 
trum at wavelengths of 127.5 nm, 130.9 nm and 186.0 nm 
respectively. While the light from the quasar may suffer sig- 
nificant absorption below the Lyman limit [corresponding 
to z <1 (< 245.0 nm)] which will influence the study of the 
Ly a lines, the analysis presented in this paper is equally ap- 
plicable to the ubiquitous metal lines that are also seen in 
quasa r spectra and can arise at lower redshifts llRauch et alJ 
12002ft . 



3 RESULTS 

3.1 Line Strength Maps 

If the absorbing material was uniform over the small scales 
probed by microlensing, then the flux in the absorption line 
and the flux in the continuum would rise and fall together 
as the source moved across the source plane. In this case, 
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Figure 6. Physical size occupied by the rays which land in the 
main region of the source plane, as a function of redshift (between 
the lens and the source, 0.0394 < z < 1.695). This relationship 
determines what area of the cloud plane is relevant for the rays 
which eventually land in the main square in the source plane. 
This explains why structure in the cloud plane is "blown up" in 
the line strength map more for clouds closer to the lens. 



the line strength map would be uniform, and a time series 
of measured values of the absorption line strength relative 
to the continuum would also be flat. Also, the distribution 
of magnifications [probability density function (PDF) for 
the magnification of a source placed at a random position 
in the source plane, see IWambseanssl dl992l ); lLewis fc Irwinl 
(1995)] for the flux in the line would be the same as that 
for the continuum, but dimmed by a constant absorption 
value. Figures [7| [S]and [5] show line strength maps for all 
of the cases investigated in this paper, together with mag- 
nification distributions for the continuum and the flux in 
the absorption line. The magnification distributions for the 
absorbed part of the spectrum are all essentially the same, 
with the most likely value having the flux in the absorption 
line about 1 magnitude fainter than the flux in the contin- 
uum. The discrepancies in the faint end for the larger cloud 
sizes are likely to be due to sampling error, with only a small 
number of clouds affecting the final magnification map. 

For the A = 0.2 and A = 0.5 cases, the structure in the 
line strength map is complex. The caustic network from the 
microlensing is evident, so rapid brightness changes due to 
the source crossing a caustic will usually be accompanied by 
sudden changes in line strength, albeit not as drastic as the 
sudden changes in the continuum brightness. Since the rays 
which landed in the relevant area of the source plane occupy 
a smaller region of space for lower A, only a small region 
of the cloud plane gets imprinted on the line strength map. 
Therefore, for a particular cloud size, faster variations in the 
line strength occur for clouds close to the source. The effect 
of the external shear in stretching both the magnification 
map and the cloud related structures in the line strength 
map in the vertical direction is clearly visible. 

When the absorption clouds are at A = 0.95 (close to the 
source) , the imprint of the caustic structure becomes much 
less obvious in the line strength map. The line strength map 
in this case approximately resembles the cloud structure, so 
if A is sufficiently high, the absorption can be approximated 
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Figure 7. Line strength maps and magnification PDFs for the 
"small" cloud size 0.005 pc, for clouds at three different redshifts. 
Changes occur on scales much smaller than an Einstein Radius, 
with the fastest variations occurring for absorption clouds close 
to the lens. The units for the magnification PDFs are magnitudes, 
with being the mean continuum magnification. 
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Figure 8. As for Figure l7l but for the "medium" cloud size 0.02 
pc. 



as being at the source plane, and line strength variations 
directly probe the structure absorbing material. As Figure|S| 
demonstrates, this approximation is applicable over a wide 
range of redshifts. 

For the low values of A, the rays which land in the region 
of interest of the source plane have only crossed a small area 
of the cloud plane, so the features in the line strength map 
are quite large. In contrast, the A = 0.95 line strength map 
shows that fast variations would occur in this case. 

For the large 0.1 pc clouds (Figure[5J, the line strength 
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Figure 9. Line strength maps and magnification PDFs for the 
"large" cloud size 0.1 pc. 

maps contain features that are typically several Einstein 
radii across, corresponding to time scales of several decades 
for any changes to be observed. Hence, this cloud size is ap- 
proaching an upper limit of the scales that could be probed 
by this method. The caustic structure is still present, but the 
same structure is also present in the line strength maps for 
the other cloud sizes, so this cannot be used to distinguish 
different cloud sizes. 

The fact that the same caustic network is present in the 
magnification map for the continuum and the line strength 
map suggests that a correlation may be present. For exam- 
ple, areas where the line strength would be strong may be 
associated with areas where the quasar has a high magni- 
fication, or vice versa. To investigate this, histograms were 
produced, plotting the line strength versus continuum mag- 
nification (one point for each pixel in the maps). These are 
shown in Figure 1101 An observation of a single light curve 
would sample regions from one of these distributions, pos- 
sibly making different cloud sizes distinguishable. In each 
case the magnification probability distribution for the unob- 
scured quasar source is the same. Clearly, Figure [T01 reveals 
that there is no strong correlation between the continuum 
magnification and the line strength relative to the contin- 
uum. Even though the imprint of the caustic network on 
the line strength maps also implies that rapid changes in the 
quasar brightness will be accompanied by significant changes 
in the line strength, there is no tendency for the changes to 
be correlated. Therefore, while spectral monitoring of the 
quasar during a high magnification event may make line 
strength variations measurable, they do not provide much 
information on the nature of the absorption clouds. Longer 
term spectral monitoring would be necessary to achieve this. 

3.2 Time Scales 

The usual way to quantify the time scale for microlens- 
ing variations is by computing t he autocorrelation func- 
tions of a sample of light curves JSeitz fc Schneiderlll994t 
ISeitz. Wambseanss. fc Schneider! Il994h . In this case, it is 
the line strength variations which are being investigated, 



Table 1. Time scales for line strength variation as a function 
of cloud size and cloud position. The units are Einstein Radii, 
which take approximately 9.5 years to cross. 





A = 0.2 


A = 0.5 


A = 0.95 


Small 


0.33 


0.24 


0.09 


Medium 


1.14 


0.81 


0.37 


Large 


3.9 


6.4 


1.7 



so the line strength map was used instead of the magnifica- 
tion map. We used a slightly different approach to measure 
the typical time scale for line strength variations as a func- 
tion of A and the size of the absorbing clouds. For a given 
line strength map, a sample of pairs of points a fixed dis- 
tance apart was selected at random, and the mean difference 
< Li — 1/2 > in the line strength values between the two 
points was calculated. 

For example, if the pairs of points are close together 
relative to the size of the variations, the line strength is 
similar for those two points, so the mean difference is small. 
In this way, a plot of mean difference versus separation was 
produced. The separation at which the mean difference is 
half of its maximum value was taken to be the measure of 
the length scale of features in the line strength map. 

Using the above method, the time scales for the varia- 
tions for each case presented in this paper were calculated 
(see Table [y. The units in the table are Einstein Radii; for 
the assumed transverse velocity of 600 km s _1 for the lens- 
ing galaxy, an Einstein Radius is crossed in about 9.5 years 
(see Section \'2. It . 

An examination of Table reveals two clear trends. 
Firstly, as expected, variability in the line strength is most 
rapid for the smallest clouds. Additionally, it is apparent 
that the time scale of variability for a fixed physical cloud 
size decreases with increasing A, and hence redshift. This 
implies that for a single microlensed quasar possessing iden- 
tical structure throughout the Ly a forest, the line variability 
will be most rapid for those absorption lines resulting from 
clouds closest to the quasar (remembering the non-linear 
relationship between A and z presented in Figure 01. 

The fastest variations occur for small clouds and clouds 
close to the source, with the variations being observable over 
a small number of decades for both the small and medium 
cloud sizes at any position between the lens and the source. 
The anomalous value for the largest clouds is probably due 
to a sampling effect, as the source plane was not large enough 
to be affected by many clouds. 

For this effect to be measurable, there must be a signif- 
icant change in the line strength over a time scale which is 
not too long. An absolute change of 0.25 in the line strength 
ought to be detectable even in low quality spectra. Recall 
that our measure of line strength is the ratio of the flux at the 
bottom of the absorption line to the flux in the nearby con- 
tinuum. We measured the distribution of the waiting times 
for a change of 0.25 in the line strength. This was done by 
selecting a random starting point in the line strength map, 
and heading off in a random direction, moving until the line 
strength has changed by 0.25. The distribution of waiting 
times for this to occur is shown in Figure ITT1 with the corre- 
sponding medians in Table|5| As expected, the waiting time 
tends to be shortest for small clouds close to the source, and 
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Figure 10. Two dimensional histograms of the line flux relative to the continuum (y axis) versus the continuum magnification (x axis) 
for each of the cases presented in this paper. The greyscale indicates the log of the relative probability of observing a particular continuum 
magnification and line strength; a strong correlation would appear as a narrow white ridge line. There are no strong correlations between 
the line strength and the conti nuum magnification, despite the sa me pa tterns being apparent in the maps. The dip in the magnification 
distribution for the continuum iLewis &; Irwii]lll995t iRauch et al.|[l993) is evident in these plots. 



Table 2. Medians of the waiting time distributions in Fig- 
ure llll in units of Einstein Radii (about 9.5 years). Note that 
the results for the large clouds are underestimated because the 
high tails of the distributions are were neglected. The peak of 
the distributions are evident though. 





A = 0.2 


A = 0.5 


A = 0.95 


Small 


0.86 


0.76 


0.26 


Medium 


1.80 


1.60 


0.70 


Large 


3.07 


3.27 


2.07 



the effect is observable in a reasonable time scale for all but 
the largest clouds. 

In an unlensed quasar, the transverse velocities of the 
quasar and the Ly Q clouds would also produce variations in 
the strength of the absorption lines. For a transverse velocity 
of 600 km s~ , the time taken to traverse 1 pc of the cloud 
plane (with clouds at redshifts z = 0.049,0.077 and 0.53) 
are 1700, 1750 and 2500 years respectively. Clearly, these 
are much longer than the corresponding time scales for the 
lensed system. For example, even for the small clouds of size 
~ 0.005 pc, the time scale is > 8.5 years, approximately 
a factor of 10 slower than the microlensed case. Thus, mi- 
crolensing provides the opportunity to study the small scale 
structure on a reasonable time scale. For the large clouds, 
the effect is even more drastic. The crossing time of a cloud 
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Waiting Time for 25% Change (Einslein Radii) 

Figure 11. Distribution of the waiting time for a 25 per cent 
change in the line strength. The x-axes are all on the same scale, 
with the time units being Einstein Radii or about 9.5 years. The 
medians of these distributions are shown in Table 121 



is ~ 2000 years, roughly 3 orders of magnitude longer than 
in the microlensed case. 
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4 CONCLUSION 

The cosmological distribution of neutral hydrogen gas leaves 
an imprint on the spectrum of background quasars, in the 
form of a series of absorption lines. This phenomenon is 
thought to be related to the large scale flow of gas and for- 
mation of galaxies and clusters. While strong constraints 
have been placed on the larger scale properties of the Ly a 
clouds, virtually nothing is known on smaller scales, where 
gas may possess similar turbulent properties to the local in- 
terstellar medium, or be collapsing to form stars. In this 
paper, we have examined the influence of gravitational mi- 
crolensing on the paths taken by light rays impinging on a 
distribution of Ly a clouds lying between the lens and the 
source. For a single macro-image, the population of micro- 
images samples differing absorption regions, with the net 
absorption seen (in a single line) being the average of that 
experienced by each ray. As the stars move, the configura- 
tion of micro-images changes, sampling different regions of 
the absorption clouds and causing temporal variability of 
the absorption line strength. 

In this study we considered clouds of three different 
sizes and at three different locations between the lens and 
the source. Quantitatively, it was shown that the fastest vari- 
ations in the absorption line strength occur for the smallest 
clouds, and for clouds close to the source. We also found 
that, due to the geometry of the universe, the variations 
in the absorption lines in the spectrum of Q2237+0305 due 
to clouds at redshifts > 0.5 would have approximately the 
same effect as identical clouds which are just in front of the 
source. At all redshifts, the caustic structure evident in the 
microlensing magnification map is also present in the ab- 
sorption line strength map. These discontinuities can result 
in very rapid changes in the absorption line strength, with 
a corresponding change in the continuum brightness. How- 
ever, a more detailed examination of the of the relationship 
between the continuum magnification and the line strength 
shows no apparent correlation between the two. 

An analysis of the variability timescales reveals that sig- 
nificant changes in line strength due to the effects of gravita- 
tional microlensing would be visible on time scales ranging 
from years to decades, dependent on the size on the size 
of the absorption clouds and their location behind the lens; 
this is at least an order of magnitude faster than clouds 
simply drifting across the line of sight to a distant quasar. 
Hence, in order to measure this effect, long term spectro- 
scopic monitoring of the individual images of microlensed 
quasars is required. Such observations have already been 
proposed for studies of quasar structu re jAb aias et al J2002t 
IShalvapin et al]l2002t iLewis k, Ibatall2004l) . 

This study has only considered Ly a clouds in image 
A of Q2237+0305, which is a positive parity image. Other 
images of Q2237+0305 have negative parity, corresponding 
to one of the curves in Figure |S] passing through zero. Fur- 
thermore, this study has assumed that the absorbing clouds 
possess no kinematic substructure, and hence the resultant 
variability applies to the total strength of the line; for this a 
simple equivalent width measure in relatively low resulution 
spectra would suffice. Nevertheless, any real distribution of 
gas will exhibit some kinematic substructure in the form of 
rotation, infall, or simple clumping. This would result in not 
only variability in the total strength of an absorption line, 



but also in its profile [see jLewis fc Bellelll99SM . To detect 
line profile variability, higher resolution spectra would be re- 
quired. We will discuss these issues in future contributions. 
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